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ABSTRACT 

We have carried out an Hi stacking analysis of a volume-limited sample of ~ 5000 
galaxies with imaging and spectroscopic data from GALEX and the Sloan Digital 
Sky Survey, which lie within the current footprint of the Arecibo Legacy Fast ALFA 
(ALFALFA) Survey. Our galaxies are selected to have stellar masses greater than 
1O 1O M0 and redshifts in the range 0.025 < z < 0.05. We extract a sub-sample of 
1833 "early-type" galaxies with inclinations less than 70°, with concentration indices 
C > 2.6 and with light profiles that are well fit by a De Vaucouleurs model. We then 
stack HI line spectra extracted from the ALFALFA data cubes at the 3-D positions of 
the galaxies from these two samples in bins of stellar mass, stellar mass surface density, 
central velocity dispersion, and NUV-r colour. We use the stacked spectra to estimate 
the average Hi gas fractions Mhi/M* of the galaxies in each bin. Our main result is 
that the Hi content of a galaxy is not influenced by its bulge. The average Hi gas 
fractions of galaxies in both our samples correlate most strongly with NUV-r colour 
and with stellar surface density. The relation between average Hi fraction and these 
two parameters is independent of concentration index C. We have tested whether the 
average Hi gas content of bulge-dominated galaxies on the red sequence, differs from 
that of late-type galaxies on the red sequence. We find no evidence that galaxies with 
a significant bulge component are less efficient at turning their available gas reservoirs 
into stars. This result is in contradiction with the "morphological quenching" scenario 
proposed by Martig et al. (2009). 
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1 INTRODUCTION 

Galaxies have long been know n to follow well-ordered se- 
quences in many properties fsee lRoberts fc Havnesfl994l . for 
a review). In the simplest classification scheme, galaxies di- 
vide into spirals and ellipticals. In more complex schemes, 
these classes are partitioned into additional sub-classes. Ir- 
respective of the exact classification scheme, there are clear 
systematic trends in the bulge-to-disk ratios, surface bright- 
nesses and in the concentration of light, all of which increase 
from spirals to ellipticals. Star formation rates and gas con- 
tent decrease along this same sequence. 

With the advent of large spectroscopic surveys of nearby 
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galaxies, such a s the Sloan Digital Sky Survey (SDSS, 
I York et al.l |2000T ). the relationships between galaxy prop- 
erties that can be derived from the combination of opti- 
cal imaging and spectroscopy, e.g. stellar mass, size, con- 
centration index, star formation rate, metallicity, dust con- 
tent, have now b een systematized and quantified in consider- 
able detail (e.g., Kauffmann et al. 2003bl; Brinchmann et all 



12004 iTremonti et all 12004 iBaldrv et al.ll2004 120061 ) . Or 

understanding of how the neutral gas content of galaxies 
relates to other galaxy properties lags far behind. 

The cold gas content of a galaxy is known to vary 
strongly with colour and star formation rate. The connec- 
tion between gas content and galaxy morphological type re- 
mains unclear. Whereas star-forming spiral galaxies almost 
always contain Hi gas, the Hi content of early-type galaxies 
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is considerably more difficult to predict. Some Es and SOs 
have neutral atomic hydrogen gas content similar to those 
of Sb-Sc type spirals , while others contain s everal orders 
of magnitude less Hi (jRoberts fc HavnesHl994r ). It has been 
speculated that these variations may be an in dication that 
the H i gas in ellipticals has an external origin ijKnapp et al.l 
Il985l ). but direct proof of this conjecture is still lacking. 

Because the Hi is on average more difficult to detect in 
early-type galaxies, the samples discussed in the literature 
have generally been quite small. Some of the largest system- 
atic studies o f gas in early- type galaxies were carried out 
in the 1980's. Knapp et al.l l|1985l ) analyzed a sample of 152 
nearby elliptical galaxies, of which 23 were detected in the 
Hi line. These authors studied the distribution of the quan- 
tity Mhi/Ls for these systems. In contrast to spiral galax- 
ies, where the distribution N(Mhi/Lb) has a well-defined 
mean value and a small dispersion, Mm /Lb spans a wide 
range in ellipticals. IWardle fc Kna pp ( 1986) extended this 
work to SO galaxies and found that the Hi detection rate 
was twice as high com pared to ellipticals. Some years later, 
iBreeman etUI l| 19921 ) carried out a study of the interstellar 
components of 467 early-type galaxies in the Revised Shap- 
ley Ames Catalogue and again reported a trend of increasing 
neutral gas content from E to Sa. These authors suggested 
that the cold gas in early-types is associated mainly with 
disks and not with the bulge components of these galax- 
ies. However, more recent studies that have ma pped Hi in 
nearby early- type galaxies (|Morganti et alj|2006l ) have con- 
cluded that the Hi can be organized in a variety of different 
configurations, e.g. in regular disks, in clouds, in rings, or 
even in tidal tail-like structures. 

One major problem that has plagued our understand- 
ing of cold gas in early-type galaxies is that the available 
Hi data have been inhomogeneous. Large area, blind Hi sur- 
veys suc h as the Arecibo Leg acy Fast ALFA Survey (AL- 
FALFA; [GjOTanenTeT^ Hool) offer uniform coverage over 
large regions of the sky and allow one to construct com- 
plete, unbiased samples of Hi-selected galaxies. However, 
these surveys are shallow, and do not in general detect gas- 
poor early-type galaxi es. This limitation has been po inted 
out in recent paper s bv ldi Sereeo Alighieri et aD (j2007r ) and 
iGrossi et"akl |2009l 'l. which used ALFALFA data to study an 
unbiased sample of early-type galaxies in the Virgo cluster 
region. They were able to compare the Hi content of early- 
type galaxies drawn from field and group environments, but 
their average detection rates were much smaller than earlier 
studies based on incomplete and inhomogeneous data. 

The Hubble morphological classification scheme is 
based on the optical appearance of a galaxy which in turn 
depends both on structural properties such as its bulge-to- 
disk ratio, and on star formation rate. If we wish to un- 
derstand the physical processes that regulate the neutral 
gas content of galaxies, it is preferable to analyze the ef- 
fects o f star formation a nd galaxy structure separately. Re- 
cently, iHelmboldd (|2007n studied 30 E and SO galaxies with 
signs of recent or ongoing star formation and concluded that 
such systems are more gas rich than E/SO galaxies with old 
stellar populations. The availability of sizes, surface bright- 
nesses and parameters measuring the concentration of the 
light for samples of millions of galaxies made possible due to 
recent advances in large scale CCD surveys and automatic 
image processing techniques, enable a new approach to un- 



derstanding the interplay between stars and gas in early- 
type systems. 

In this paper, we make use of stacking techniques to 
analyze whether the average Hi gas fraction of a galaxy is 
affected by the presence of a significant bulge component. 
Stacking has now become a common tool to constrain the 
statistical properties of a population of objects that lack in- 
dividual detections in a survey; by co-adding the signal from 
many objects with known sky positions and redshifts, the 
background noise can be decreased and one can recover the 
average Hi flux of the ensemble. Stacking techniques have 
been applied to a wide variety of different astro physical data. 
Exam ples include studies of faint radio AGN (jHodee et alJ 
2009), studies o f star formation in high redshift Lyman 
Break Galaxies (jCarilli et al.l 120081 ) , studies of the intra- 
clust er light using stacked optical images (e.g. JZibetti et aD 
120051). and the sof t X-ray properties of high redshift quasars 
i Shen et alJ I2006T ). Stacking has been applied to Hi data 
as well, with the purpose of studying the Hi properties of 
gas at redshifts that are currently not well probed by ex- 
isting radio telescopes (individual detections of Hi emission 
reach redshifts z^0.25 and require extremely long integra- 
tions, stacked 
non-detections in different regions of a z=0.06 cluster to in- 
vestigate en vironmental effec t s. Mo re recently, Hi stacking 
was used by IVerheiien et al.l (12007 ) to probe the Butch er- 
Oelmer effect at z~0.2, and bv lLah et aD (|2007l . l2009h to 
attempt to constrain the Hi content of star-forming galaxies 
at z~0.24 and galaxies around a cluster at z~0.37. 

Here, we use ALFALFA survey data to constrain the 
average Hi gas fractions of an unbiased sample of massive 
early-type galaxies. We study how the Hi content depends 
on parameters such as stellar mass, stellar mass surface 
density, concentration index, central velocity dispersion and 
UV/optical colour. The paper is structured as follows. In 
section [2] we describe the samples considered in this pa- 
per. The stacking analysis is described in section [3] First 
we study Hi gas fraction scaling relations for a complete 
sample of galaxies and then we compare our results with 
previous work (section [4}. The analysis of the Hi properties 
of early- type galaxies is presented in section [5] All distance- 
dependent quantities in this work are computed assuming 
Q. m =0.3, A =0.7 and H =70 km s^Mpc" 1 . 



2 SAMPLE SELECTION 

Measured redshifts are essential if we wish to recover an 
accurate estimate of the mean Hi content of a popula- 
tion of galaxies using stacking techniques. Our sample is 
drawn from th e "parent sample" of the GALEX Arecibo 
SDSS Survey ijCatinella et al.l l2010l . hereafter GASS-1), 
which is a sample of 12006 galaxies with stellar masses 
greater than 10 10 AfQ and redshifts in the range 0.025 < 
z < 0.05 selected from the SDSS main spectroscopic sam- 
ple. The parent sample galaxies are located in the in- 
tersection o f the footprints of the Data R elease 6 of the 
SDSS (DR6; lAdelman- McCarthy et al.ll2008h. the projected 
GALEX Medium Imaging Survey fMIS; lMartin et al.ll2005h 
and ALFALFA. The average uncertainty in the SDSS spec- 
troscopic redshifts is extremely small (0.0002). As discussed 
in GASS-1 (see section 3 for further details), the stel- 
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Figure 1. Distributions of (a) stellar mass, (b) redshift and (c) NUV-r colour (corrected for Galactic extinction only) for galaxies in 
sample A. The black solid histograms represent the whole sample A, while the red, dashed histograms show the distributions for the 
sub-sample of galaxies with ALFALFA detections. Panel (d) shows the colour-magnitude diagram of sample A galaxies (gray dots) and 
the sub-sample with ALFALFA detections (red dots). 



lar masses are d erived from SDSS photometry following 
ISalim et al.l l|2007l ) and have typical errors smaller than 30%. 
The choice of redshift range is determined by sensitivity lim- 
its in the Hi observations and by the need to avoid redshift 
ranges where radio frequency interference (RFI) is a prob- 
lem. 

The acquisition of ALFALFA data is on-going. Part 
of the data have already been catalogued and are avail- 
able to the publ i c (for t he ALFALFA sky c o vered by SPSS: 
Giovanelli et alj 120071 : iKent et alj |200S| ; IStierwalt et~aT] 
20091 ). The ALFALFA 40% dataset to be released in late 
2010 (Martin et al. in preparation; Giovanelli et ai, in prepa- 
ration) includes the following SDSS sky regions: 7.5 h < 
Q2000 < 16.5 h, +4° < 5 2 ooo < +16° and +24° < 5 20 oo < 
+28°, and 22 h < a 2 ooo < 3h, + 14° < 5 20 oo < +16°. 

Within the same sky region, the GASS sample contains 
5350 galaxies; these constitute the parent sample for our 
study. We discard some objects because they have poor qual- 
ity Hi data (see section OS), so that the final GASS-SDSS- 
ALFALFA sample (which we call sample A) is composed of 
4726 objects. Of these, 23% are galaxies with reliable AL- 
FALFA detections (i.e. objects corresponding to ALFALFA 
detection codes 1 or 33). Panels (a), (b) and (c) of Fig- 

1 As discussed in lGiovanelli et al.1 i2005T) , ALFALFA HI line de- 



ure[T]show the stellar mass, redshift and NUV-r colour dis- 
tributions of the galaxies in sample A. The NUV-r colours 
have been corrected for foreground Galactic extinction, but 
not for internal extinction. The ALFALFA detection rate is 
close to 23% for each stellar mass bin in Figure [T] but is 
clearly biased to blue-sequence objects (with NUV— r <3.5) 
as shown in panel (c). Finally, in panel (d) we plot NUV-r 
colour versus absolute r-band magnitude M r for the sample, 
with black dots representing the full sample and red points 
the galaxies detected by ALFALFA. Once again we see that 
the galaxies detected by ALFALFA are almost exclusively 
found on the blue sequence. 

2.1 Galaxy Parameters 

The optical parameters we use are derived from the MPA- 
JHU SDSS DR7 release of spectrum measurements or from 
Structured Query Language (SQL) queries to the SDSS DR7 

tections are coded into two categories: Code 1 detections have 
a peak signal-to-noise ratio greater than 6.5 and are reliable at 
greater than 95% confidence; Code 2 detections, referred to as 
"priors", have a lower signal-to-noise ratio between 4.5 and 6.5 
but an optical counterpart at the same known redshift. Their re- 
liability is estimated to be greater than 85%. 
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Figure 2. As in Figure 1, except for the ETG sample. 



database servei0- The UV para meters are ex t racted from 
the GALEX UV photometry by IWang et all (|2009l ). The 
reader is referred to section 5 of GASS-1 for more detailed 
descriptions. 

The parameters used in this paper are the following: 1) 
stellar mass M*, 2) stellar mass surface density fj,*, defined 
as ^i*=M*/(27ri?5 0jZ ), where -Rso.z is the Petrosian radius 
containing 50% of the flux in z-band in units of kpc, 3) the 
concentration index C—Rgo/Rso, where Rqo and R50 are 
the radii enclosing 90% and 50% of the r-band Petrosian 
flux, 4) and NUV-r colour. As already explained, the 
colour is corrected for Galactic extinction only. Corrections 
for internal dust-atte n uatio n are discussed and applied 
in ISchiminovich et al] |2010), where they study the star 
formation properties of the GASS sample. 

In addition to these quantities, we extract the following 
photometric parameters: 

(i) The inclination i to the line-of-sight is evaluated ac- 
cording to: cosi = b/a, where a and b are the semi-major 
and semi-minor axes from the r-band exponential fit, respec- 
tively (b/a for the exponential fit is tabulated as expAB_R). 

(ii) The likelihood parameters (Inldevjr and Inlexpjr from 
SDSS), which indicate how well the de Vaucouleurs and the 



exponential models fit the one-dimensional r-band light pro- 
file of the galaxy 

(iii) The central velocity dispersion (j 1( / g is derived from 
the SDSS parameter vdisp. These velocities are evaluated 
with the "direct fitting" methodjfl using spectra measured 
within the 3"-diameter fiber aperture. Only values between 
70 and 420 kms -1 are rel iable. We the n correc t vdisp for 
aperture effects. Following iGraves et all (|2009al ) and refer- 
ences therein, we scale the fiber velocity dispersion to that 
at 1/8 of the effective radius: ai/g = vdisp - (r/it>/|r ) ' 04 , 
where r/it, = l".5 and r is the circular galaxy radius, de- 
fined as r = R e \Jbja. R e is the r-band de Vaucouleurs 
radius (tabulated as devRjr), and a, b are now the major 
and minor axis from the de Vaucouleurs fit (b/a is tabu- 
lated as devabjr). In general, this correction is small, ~ 5%. 
(Note that for galaxies included in sample A, R e has a mean 
value of 8.2".) 



2.2 ETG sample selection 

In this paper, we have chosen to define "early-type" galaxies 
purely in terms of their structural properties, without re- 
gard to their stellar populations or star formation rates. We 
note that our definition is in contrast to some definitions of 



2 See http://www.mpa-garching.mpg.de/SDSS/DR7/ 
http: / /cas.sdss.org./dr7/en/tools/search/sql.asp 



and 



3 See http://www.sdss.org/DR6/algorithms/veldisp.html for 
more discussion 
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"early-types" in the liter ature, which have exc l uded galaxies 
with emission lines fe.g., lBernardi et alJliool lG raves et al.l 
l2009al ). Our goal will be to explore the extent to which the 
presence or absence of a significant bulge component influ- 
ences the Hi content of a galaxy, so we do not wish to bias 
our conclusions by selecting against bulge-dominated galax- 
ies with emission lines. 

Starting from sample A, we extracted a subset of early- 
type objects {ETG sample) with the following properties: 

• concentration index C—R30/R50 > 2.6; 

• the likelihood that the light profile is fitted by a de 
Vaucouleurs model is greater than it is by an exponential 
one; 

• inclination less than 70° (this cut rejects the most in- 
clined systems). 

We note that the C parameter has been shown to be 
an excellent indicator of the bulge-to-total ratio (B/T) de- 
rived from_foll2j ; dimension bulge/disk decomposition anal- 
ysis (|Gadotti| [2009). A cut at C>2.6 res tricts the sample 
to ga laxies with B/T > 0.4 (see Fig.l of IWeinmann et al.l 
120091 ). For the present work we choose our default cut at 
C > 2.6, but we also experiment with cuts at larger values 
of C. 

Application of our default cut leads to a final ETG sam- 
ple consisting of 1833 objects. The properties of this sample 
are shown in Figure [2] The solid histograms represent the 
entire sample; the red dashed histograms the sub-sample 
detected by ALFALFA. The stellar mass, redshift and the 
NUV-r colour distributions are shown in panels (a), (b) and 
(c), respectively. The average ALFALFA detection rate for 
early-type objects is smaller (~ 9%) than for sample A. The 
colour- magnitude diagram is shown in panel (d). Most of the 
early-type targets lie on a well-defined red sequence. Some 
objects do scatter bluewards of the red sequence. Such ob- 
j ects may be star-form ing, transitional or Seyfert galaxies 
l|Schawinski et al]|2007l ). By selecting targets based on con- 
centration index and inclination, we may also include objects 
with some disk component. In Figure Owe presents lxl ar- 
cminute SDSS postage stamps of a randomly selected set of 
galaxies from our ETG sample, that have been ordered by 
increasing C. 



3 ALFALFA DATA STACKING 

The on-going ALFALFA survey is scanning 7000 deg 2 
of the high galactic latitude sky over the veloc- 
ity interval «[kms _1 ] ~ [-2500; 18000] (i.e. out to 
z~0.06) u sing a simple "m i nimu m-intrusion" drift scanning 
technique (|Giovanelli et al] I2005T ) that exploits the seven- 
horn Arecibo L-band feed array (ALFA). For each of the 
seven beams, spectra are recorded separately for the two 
orthogonal polarizations, providing two independent sam- 
ples. The 2-D drift scan data from both polarizations and 
all beams covering a given portion of the sky are combined to 
form 3-D cubes of dimension 2.4° x 2.4° on the sky and 5500 
kms -1 in velocity "depth". We refer the reader to Figure [H 
where a schematic representation of one such data-cube is 
shown. We will refer to this depiction throughout the paper. 

The raw spectral resolution of the ALFALFA data, be- 
fore smoothing, is ~5.5 kms -1 and the angular resolution 



is ~3.3'x3.8' (corresponding to the FWHM of each ALFA 
beam). The data-cubes or "grids" are constructed from the 
drift scan data so that each spatial pixel is 1' on a side. The 

2- D drift scans are flux calibrated using the real-time noise 
diode injection scheme; the position and flux scales of the 
final grids are updated us ing fits to the ma ny radio contin- 
uum sources they co ntain dKent et al j200Sfl . The ALFALFA 
processing scheme l|Giovanelli et aLl2005f ) retains all indi- 
vidually recorded spectra for each ALFA beam and each of 
its polarizations separately; no filtering is performed to dis- 
card bad channels. In order to deal with poor quality data, 
the 2-D data for each beam and polarization are visually 
inspected to flag bad records and frequency channels. The 

3- D grid construction then proceeds with knowledge of the 
flagged 2-D pixels. Data quality in the final grid then may 
be limited by the absence of data in channels flagged as con- 
taminated by RFI, weighted as poor by occasional instru- 
mental problems or missing entirely because of incomplete 
sky coverage. In order to account for these effects, each 3-D 
pixel of the data-cube is assigned a quality weight w during 
the data reduction process, which is a number ranging from 
(unusable data) to 20 (good data), computed according 
to the availability and quality of data contributing to each 
pixel. The retention of the full dataset and the construction 
of the accompany "weights map" allows us to judge whether 
or not adequate data exist for each target so that it can be 
meaningfully included in the stacking process. 



3.1 Creating a catalogue of Hi spectra 

We extract a spectrum for each galaxy in the sample. 
As discussed, all our targets are selected from the SDSS 
spectroscopic survey, so we know both their positions on 
the sky and their redshifts. We first select the ALFALFA 
data-cube which contains the target and then follow a 
procedure that includes the following steps: a) spectrum 
extraction; b) rms evaluation, c) final quality check. 

a. Spectrum extraction 

The signal from each target is integrated over a region of 
the data cube centered on its 3-D position. Because noise 
increases with the square root of the integration area, in- 
tegrating over too large a region lowers the quality of the 
spectrum without increasing the signal. Our GASS targets 
are always smaller then the ALFA beam (the mean i?go for 
sample A is 10"), so we simply integrate over a sky region of 
4' x 4' . In Figure [4] we illustrate how we extract spectra at 
two different positions in the sky inside the same data-cube. 
The coloured regions indicate where the spectra would be 
evaluated. 

The Hi spectrum is a histogram of flux density S as a 
function of velocity. For each velocity channel v, the cor- 
responding flux density S v is obtained by integrating the 
signal s v (x,y) over the spatial pixels centered at the target 
galaxy position, as observed by a radio telescope of beam 
response pattern B: 



S v [mjy] 



Yl x YlyS v {x, y) 



where x, y are the sky coordinates (the two polarizations 
are kept separated). The expression above means that the 
spatially integrated profile is obtained by summing the sig- 
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Figure 4. A schematic representation of the fully processed ALFALFA 3-D data-cube. The data cubes are 2.4° X 2.4° in size and about 
5500 kms -1 in velocity range (25 MHz in frequency). The raw spectral resolution is ~5.5 kms -1 ; the angular resolution is <~4'. For each 
pixel, which is a point in RA, Dec and velocity, a value of flux density is recorded. For each target in sample A, we extract a spectrum 
at a given position of the sky, over the velocity range of the data-cube which contains the source. Two examples of extracted spectra are 
shown on the right, illustrating an Hi detection (green, bottom) and an Hi non-detection (red, top). 

rials over all the spatial pixels of interest and dividing by ALFALFA beam pattern can be approximated by: 
the sum of the normalized bea m B(x,y) over the sam e pix- 
els (for a detailed discussion see lShostak fc A llen 1980). The 
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with a x = (2V21n2)~ 1 x 3.3', and a y = (2V21n2)~ 1 x 3.8' 
l|Giovanelli et alj|2005l ). 

We note that we discard any spectrum if more than 
40% of the pixels have a a quality weight w less than 
10. We also keep track of the three strongest continuum 
sources in an area covering 40'x40' around each source; 
strong continuum sources can affect our spectra by creating 
standing wavefl 

b. Rms evaluation 

For each spectrum we need to measure the root mean square 
(rms) noise, which will later be used as a weighting factor 
when we stack spectra. The rms has to be evaluated in re- 
gions of the spectrum where there is no emission from the 
target galaxy, and we also have to avoid spectral regions 
where there are any spurious signals (e.g. residual RFI that 
we failed to flag, or Hi emission from companion galaxies). In 
order to define the spectral region that might contain galaxy 
emission, we estimate its expected Hi width as follows. 
The expected width of the Hi spectrum will depend on 
the rotational velocity of the galaxy as observed along the 
line-of-sight. We estimate the expected velocity wtf-,o for 
each target, using the the Tully-Fisher relation. Following 
iGiovanelli et~ai1 (|l997bh . we use the SDSS i-band magnitude 
(fc-corrected and corrected for Gal actic and internal extin c- 
tion as in equations 11 and 12 in IGiovanelli et ail (|l997aT )i 
to estimate wtf , and the measured inclination of the galaxy 
to derive wtf-.o- We are aware that the Tully-Fisher relation 
does not hold for all morphological types and environments. 
We do not think this is a major issue, because these veloc- 
ities are only used to estimate the region of the spectrum 
that should contain significant signal from the galaxy. 

We then fit a first order polynomial to the baseline af- 
ter excluding the region of the spectrum containing signal 
from the galaxy. This step allows us to eliminate possible 
gradients in the background. (The top right panel in Fig- 
ure [4] shows an example of a spectrum where the baseline is 
tilted). We perform a robust polynomial fit over the regions 
of the spectrum with high values of the quality factor w and 
then evaluate the rms about the fit over the same region of 
the spectrum. 

The average rms for the whole sample (for each 
polarization) is 3.6 ± 0.5 mjy After averaging the two 
polarizations, an rms of ~2.5 mjy is obtained; this is 
comparable to the average rms of 2.2 mjy evaluated for 
published, reduced ALFALFA spectra. 

c. Final quality check 

After we have extracted the spectra, we visually inspect each 
of them. We check the extraction process, and we discard 
spectra with bad baselines caused by continuum sources and 
those with possible spurious signals close to the galaxy (e.g., 
if one polarization has a significantly stronger signal than the 
other, or if there is a strong signal close to the object, which 
may arise from a companion galaxy). These cuts eliminate 
624 objects in sample A (11% of the initial sample). 



4 Standing waves are periodic fluctuations in the background 
which occur when radiation from a strong continuum source is 
multiply reflected and scattered by the telescope structure before 
reaching the receiver. 
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Figure 5. Dependence of the rms of the stacked spectra as a 
function of the number of objects co-added. The dashed line is the 
expected 1 / \/~N dependence. The flattening of the relation around 
N ~ 300 arises, as expected, where the non-Gaussian noise be- 
comes dominant. 

3.2 The Stacking Method 

We want to co-add the signals from N different sources lo- 
cated at different redshifts. First we shift each spectrum to 
the galaxy rest frequency, so each spectrum is centered at 
zero velocity. We stack together the spectra Si (i=l,..N) us- 
ing their rms as a weight, a standard approach in stacking 
analysis. In doing this, the final spectrum S 3 tack would be: 

Jatack — x^N K - 1 ) 

^i=0 Wi 



The stacking of the spectra is done separately for each polar- 
ization. Note that ALFALFA is a blind survey which scans 
the sky uniformly, so the rms for most spectra are similar. 
If the noise of the input targets is purely Gaussian, the rms 
of the stacked spectrum is expected to decrease as l/\ // ~N, 
where N is the number of objects co-added. In reality, in 
addition to the Gaussian noise there are likely to be system- 
atic components, for example standing waves (we discarded 
spectra dominated by standing waves, but a weak residual 
signal could still remain). Because of these additional noise 
sources, we expect the rms to approach a lower limit as the 
number of co-added spectra becomes very large. 

In Figure [S] we show how the measured rms of the 
stacked spectra decreases as a function of the number of co- 
added objects. We have stacked increasing numbers of ran- 
domly selected galaxies from sample A, and for each stack 
we evaluated an rms as described above. These measure- 
ments are shown as dots, while the dashed line shows the 
expected average rms value assuming Gaussian noise, i.e. 
2.5 /\fN mjy. As expected, we see the trend flatten as N 
approaches values of around 300, where the non-Gaussian 
noise becomes dominant. The rms continues to decline as N 
increases, but at a slower rate. 

After stacking, we manually process each ensemble 
spectrum. The reduction process includes the following 
steps: i) we average the two polarizations, ii) A default Han- 
ning smoothing is performed on the spectrum. Depending 
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on the signal-to-noise, we may apply a boxcar smoothing to 
decrease the noise if the signal is marginal. Note that we 
never average over more than nine channels (corresponding 
to Av ~ 50 kms -1 ). Hi) Finally, the baseline is subtracted. 
We note that the baseline of a stacked spectrum is already 
almost flat, because the different noise features tend to can- 
cel out when averaging many spectra. 

Figure [6] shows some examples of stacked spectra, ob- 
tained by stacking sample A galaxies in five different stellar 
mass intervals. The vertical axis in each panel shows flux 
density in mjy, while the horizontal one shows velocity in 
kms -1 . Since we shifted each object to the rest-frequency 
of each galaxy, the stacked signal is centered at zero ve- 
locity. For each stacked spectrum, the mass range and the 
number of objects stacked are reported, as well as the signal- 
to-noise ratio of the detected Hi signal. As explained above, 
when the signal-to-noise is low, we smooth the spectrum 
over up to nine channels as done, for example, for the last 
spectrum in Figure [6] In the left column of the Figure, we 
show spectra obtained by stacking ALFALFA detections and 
non-detections. In the right panel we show the spectra ob- 
tained by stacking only the non-detections, demonstrating 
that the stacking process recovers a signal even if no in- 
dividual galaxy is detected (dotted lines in Fig. [6] indicate 
the edges of the signal). Notice that the rms decreases with 
increasing number of co-added objects. The width of the 
profile is smaller for less massive objects, as expected since 
they have on average lower circular velocities. 

If we recover a signal in the stacked spectrum, we mea- 
sure the integrated emission between the two edges of the 
Hi profile (see dotted lines in Figure [6]) , which are defined 
manually for each spectrum. We expect the signal to be sym- 
metric around zero velocity. Even if in some cases the SDSS 
redshift is slightly off-set with respect to the Hi emission, 
the discrepancy will be random and will cancel out when 
co-adding multiple spectra (Note that the SDSS redshift un- 
certainty is 0.0002, which corresponds to ~60 km s" ). 

W e evaluate a signal-to-noise ratio following ISaintongel 
i|2007l ). We define an Hi detection if the S/N is greater than 
6.5. If there is no detection, we evaluate an upper limit, 
assuming a 5 a signal with a width of 300 kms -1 , smoothing 
the spectrum to 150 kms -1 . (The width of 300 kms -1 is 
chosen because it corresponds to the peak of the distribution 
of velocity widths for galaxies in the GASS survey.) 



3.3 Evaluating Hi gas fractions 

Our aim in stacking spectra is to characterize the average 
Hi content of a given sample of galaxies, so we are interested 
in converting our recovered signal into an Hi mass and also 
into an average Hi gas fraction. For a single object, these 
two quan tities are well defined. The Hi mass is evaluated 
using the lRobertsl (|f963T ) formula: 

Mhi _ 2.356 x IP 5 ( Dl{z) \ 2 ( S int \ . . 

M© 1 + z V Mpc J VJykms-V' 

where Dl(z) is the luminosity distance and Si n t the inte- 
grated Hi line flux. A correction for Hi self-absorption is 
not applied to the Hi mas s, as it is likely to be negligible 
l|Havnes fc Giovanelli|[l983 , Appendix B). The Hi gas frac- 
tion is simply defined as Mhi/M*. 



For N individual detections (each with a flux measure- 
ment), the average value of their gas fractions can be defined 
by the weighted mean value: 

,M H I, f„N MHI;t \ , („N \ fn a\ 

{ nr ] = r feo -M^r" u, v / ( E< = oU, (3 - 4) 

where the Wi are defined as in equation 3.2. 

When co-adding spectra, we need to take into account 
the fact that our targets span a significant range in redshift 
(zmax = 2 • z m in), so given the same Hi mass an object 
at the lowest redshift limit contributes 4 times more signal 
than one at the upper redshift limit. Because we are stacking 
mainly non-detected spectra we do not know how much each 
galaxy contributes to the total signal. Both the mean and 
median values of redshift and stellar mass may not be repre- 
sentative. In order to weight each spectrum in a consistent 
manner, we choose to stack "gas-fraction" spectra S[, where 
the signal that we co-add is no longer the flux S [m Jy] , but 
the quantity: 

St [mjy] -> S[ [mjy Mpc 2 M^ 1 ] = Sl ' (3.5) 

The stacking is performed according to equation 13.11 With 
this approach, the stacking creates an average object with 
respect to the Hi gas fraction of the individual galaxies used 
to build the stack. In Appendix A, we discuss an alterna- 
tive approach to evaluating mean gas fractions from stacked 
spectra, and present comparisons between the two different 
methods. 

We meas ure errors on the gas frac tions using the Jack- 
knife method (|Tukevll 19771 : lEfronlll982r ), a statistical tool to 
estimate a confidence interval on a function of N measures. 
The purpose of computing jackknife errors on the mean Hi 
fraction estimated from a stacked spectrum is to ascertain 
whether or not the signal is dominated by a few outliers. 
Schematically, given a function &(xi) of N observations Xi, 
the Jackknife evaluates N partial estimates ay of the func- 
tion obtained by discarding one element per time. The jack- 
knifed estimate a* is then the average of the pseudo-values 
a* — Na — (N — Finally, the sample variance on the 

pseudo- values is: 

1 N 

{ojackf = N ( N _^ ~ «*) 2 > ( 3 - 6 ) 

This can be used to provide a confidence interval on a — a* . 



4 Hi GAS FRACTION SCALING RELATIONS 
FOR MASSIVE GALAXIES 

In this section, we characterize how average Hi gas fraction 
depends on a variety of different properties for galaxies with 
stellar masses greater than 1O 1O M0. First we study the de- 
pendence on stellar mass M*, stellar mass surface density 

concentration index C and colour NUV— r for galaxies 
in sample A. We then compare the same relations for galax- 
ies in the bulge-dominated ETG sample (§[5}. 

We begin by comparing our results with those of 
Catinella et al. (2010), who studied the correlation between 
Hi gas fraction and the properties listed above, for a com- 
plete sample of ~ 200 galaxies from the GASS survey with 
M* > 10 10 Mq. In that survey, the Hi observations reached 
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Figure 6. Examples of stacked spectra for 5 M* bins. The x-axis is velocity in kms -1 , the y-axis is Hi flux density in mjy. Since we 
shifted each object to its galaxy rest-frequency, the signal is centered at zero velocity. For each spectrum, the mass range and the number 
of objects stacked are reported, as well as its S/N ratio. Dotted lines show the boundaries of the signal, inside which we integrate the flux. 
In some spectra there are spikes/holes caused by poor quality data (note that the spectra containing the bad data were not discarded 
because the bad pixels are located away from the central regions of interest). Examples of bad regions occur in the first two rows between 
v =-1500 kms" 1 and v =-1000 kms^ 1 . Left column: stacked spectra using all galaxies. Right column: stacked spectra using only those 
galaxies that were not detected by ALFALFA. As expected, the signal is systematically lower for the right column. 
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Figure 7. Red circles: Average gas fractions of sample A galaxies derived from stacking are plotted as a function of stellar mass M, 
(top) and concentration index C (bottom). The points are plotted at the mean value of M, or C for the galaxies in the bin. The red line 
(upper panel) show a linear fit to these results, which are compared with the average Hi gas fractions from the GASS-1 paper (green 
lines). Green dashed lines show the la uncertainties on the GASS-1 estimates. Gray dots show galaxies with ALFALFA detections from 
sample A. The numbers written in the panels indicate the numbers of objects co-added in each bin (Tot), and the percentage of them 
directly detected by ALFALFA (%Det). The values plotted are reported in Table IAT1 4th and 5th columns. 
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Figure 8. Average gas fractions of sample A galaxies derived from stacking are plotted as a function of stellar mass surface density 
/i* (top) and NUV— r colour (bottom). Symbols and colours as in Figure [7] The values plotted are reported in Table |A"T1 4th and 5th 
columns. 
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a much deeper flux limit compared with ALFALFA, so av- 
erage Hi gas fractions could be derived using individual Hi 
detections, rather than by stacking. It is important to check 
whether our stacking method gives answers that are consis- 
tent with these published measurements. 

In Figures 7 and 8, we present our results for sample A. 
In each panel, gray dots show the galaxies in sample A that 
were detected by ALFALFA. Green lines show the average 
values and la confidence interval of Mhi/M* from GASS-1 
(see Figure 9 in that paper). Our own estimates of the Hi 
gas fraction derived from the stacked spectra are plotted as 
red circles, and the errorbars on these points are evaluated 
using the Jackknife method (eq. 13-fcif) . Note that the stacked 
spectra that yield the measurements plotted as red circles in 
Figure [7] (top panel) , are shown in Figure HJ] (left panel) . At 
the bottom of each panel, we record the number of galaxies 
included in each stack, as well as the percentage of galaxies 
with ALFALFA detections. In cases where the stack includes 
more than a few hundred galaxies, the errors on the mean 
Hi fraction are negligible. 

We note that we are able to measure an Hi gas fraction 
in every bin, even in the very reddest NUV-r colour bins 
where the ALFALFA detection rates are close to zero. 
The fact that the jackknife errors remain relatively small 
for these bins, indicates that the stacked spectra are not 
dominated by signal from a small fraction of the galaxies. 

We find that the average Hi gas fractions derived from 
the stacked spectra are in excellent agreement with the re- 
sults reported in GASS-1. The point-by-point agreement is 
generally within la, with the possible exception of the lowest 
stellar mass and stellar mass surface density bins, where we 
find gas fractions that are systematically lower than those 
reported in the GASS-1 paper. We caution that low mass 
galaxies were somewhat under-represented in the first GASS 
data release. Our errorbars only indicate statistical errors 
and do not account for other effects, such as cosmic vari- 
ance, so we do not think the discrepancy is significant. 

As discussed in GASS-1, the Hi gas fraction is a decreas- 
ing function of stellar mass, stellar mass surface density, 
concentration index and colour. The slopes of these rela- 
tions for our sample A are listed in the 3rd column of Table 
I All Here we only report the total change in log Mm /M* 
(AlogMm/M*) for each relationship, which allows a com- 
parison among the different properties: 



M* 
C 
M* 

NUV-r 



A log M H i /M* 
A log Mm /M, 
A log Mhi/M* 
A log Mm /M* 



= 0.20 
= 0.24 
= 0.45 
= 0.52 



The correlations between gas fraction and both M* and con- 
centration index (Figure [7|) are weak. In contrast, the aver- 
age value of Mhi/M* drops by more up to a factor of 25 when 
plotted as a function of NUV-r colour or stellar mass surface 
density. Although stacking recovers the average trends of the 
population extremely well, it does not provide any informa- 
tion on the underlying scatter. This can only be studied with 
individual gas fraction measurements, as in GASS-1. 



5 Hi STUDY OF A COMPLETE SAMPLE OF 
ETG GALAXIES 

In this section, we ask whether bulge-dominated galaxies 
with M* > 1O 1O M0 lie on the same Hi scaling relations as 
the general population of galaxies with M* > lO lo A/0. Our 
goal is to determine whether the presence of the bulge plays 
a role in regulating the rate at which gas is c onsumed into 
stars , for example by stabilizing the disk fe.g.. iMartig et ail 
2009). In order to test the role of the bulge in a clean way, 
we must account for the fact that the physical properties of 
galaxies are strongly correlated. If one selects a sub-sample 
of bulge-dominated galaxies from the parent sample A, one 
will automatically select a sample of galaxies with higher 
stellar masses, higher stellar mass surface densities, and red- 
der colours. It is therefore important to understand whether 
or not bulge-dominated galaxies differ in Hi content from 
the parent sample at fixed values of these parameters. 

Our results are presented in Figures [9] and [10] and in 
Table IA1I In all the plots, blue circles are the average gas 
fractions obtained from stacking the ETG sample and the 
errorbars are evaluated using the Jackknife method (eq. 13.61) . 
Upside-down triangles indicate the upper limit in the case 
of a stack that yields a non-detection. The red lines (circles) 
show the fits to the mean Hi gas fraction relations obtained 
for sample A. Our main result is that the average Hi gas 
fractions of bulge-dominated galaxies are significantly lower 
(by approximately a factor of 2) than those of the parent 
sample at a given value of stellar mass. A similar, but weaker 
reduction in the average Hi gas fraction is seen for the ETG 
sample when it is plotted as a function of the central velocity 
dispersion of the bulge. However, the relation between gas 
fraction and stellar mass surface density and NUV-r colour 
appears to be insensitive to the ETG cut. In GASS-1 paper, 
Catinella et al. showed that a linear combination of stellar 
mass surface density and NUV-r colour provided an excellent 
way to "predict" the Hi content of galaxies more massive 
than 10 10 M Q . Here we show that this conclusion holds true 
independent of the bulge-to-disk ratio of the galaxy. We have 
tested that this conclusion still holds if we define the ETG 
sample using more stringent cuts on concentration index 
(C > 3) or on the axis ratio of the galaxy (b/a > 0.6, which 
implies inclination lower than ~55°). These cuts reduce the 
ETG by 50% and 30%, respectively. Nevertheless, results 
shown in Figures 9 and 10 demonstrate that the average Hi 
gas fractions of these systems still lie on the same relations 
when plotted as a function of stellar mass surface density 
and NUV-r colour (cyan circles represent the further cut in 
concentration index, black ones the cut in inclination). 

In Figure ITTl we show how the average gas fractions of 
galaxies in the ETG sample vary as a function of position in 
the two-dimensional plane of colour versus stellar mass den- 
sity /i*. Bulge-dominated galaxies are mainly found on the 
red sequence, but there is a minority population with bluer 
colours. We adaptively bin the sample in two dimensions 
by recursively dividing the plane into axis-aligned rectan- 
gles. We stop dividing a region when a further split would 
lower the S/N below the detection threshold of 6.5. Figure 
[11] (bottom panel) shows the final binning used. In each bin 
the measured gas fraction (expressed as a percentage of the 
stellar mass) is reported. In Figure [TT] (top panel) we colour- 
code the (NUV— r)— fi* plane according to gas fraction. The 
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Figure 9. The dependence of the average Hi gas fraction on stellar mass M t (top) and on central velocity dispersion a (bottom) for the 
ETG sample (blue symbols). The relations found for sample A are shown in red for comparison (the fit for the M* relation, the actual 
points for a). Upside-down triangles indicate upper limits in the case of a non-detection. The numbers written in the panels indicate 
the numbers of objects co-added in each bin (Tot), and the percentage of them directly detected by ALFALFA (%Det). Gray dots show 
sample A galaxies with ALFALFA detections. We have also applied more stringent cuts to the ETG sample, as explained in the text: 
cyan circles represent a sample with C>3 (52% of the original ETG sample); black circles are for a sample with b/a>0.6 (or inclination 
lower than 55° - 68% of the original sample). The values plotted for each ETG sample are reported in Table [All 6th-llth columns. 
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Figure 10. Average Hi gas fraction dependence on stellar mass surface density fi t (top) and NUV— r colour (bottom) for the ETG 
sample. Symbols and colours are the same as described in Figure [9] The values plotted for each ETG sample are reported in Tabic |A"T1 
6th-llth columns. 
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Figure 11. Average Hi gas fraction dependence in the 2-dimensional plane of stellar mass surface density fi* and NUV-r colour for 
galaxies in the ETG sample. In the top panel the dots show individual objects, while the colors show the (interpolated) gas fractions 
measured with the stacking as a function of position in the plane (the colour scale key is included above the plot). The bottom panel 
shows the adopted binning. For each bin, the gas fraction measured from the stacked spectrum is noted. 
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Hi content decreases going from left to right (towards in- 
creasing stellar mass surface density) and from bottom to 
top (towards redder colours). The most significant variation 
is clearly along the colour direction. 

5.1 A Test of the morphological quenching 
scenario 

The idea that galaxy disks are more resistant to the forma- 
tion of bars, spiral density waves and other instabilities if 
they are embedded within a dy namically hot halo o r bulge , 
has its or igins in early work by lOstriker fc Peebles! l|l973t ). 
Recently, iMartig et al.l (|2009T ) have proposed this so-called 
"morphological quenching" mechanism as a way of explain- 
ing why present-day bulge-dominated galaxies on the red se- 
quence cease growing in stellar mass in environments where 
they continue to accrete gas. In their picture, a disk with 
similar gas content will be much less efficient at forming 
stars if it is embedded in a galaxy with a significant bulge 
component. As stated in the abstract of their paper, "our 
mechanism automatically links the colour of the galaxy with 
its morphology, and does not require gas consumption, re- 
moval or termination of the gas supply." 

To test whether the "morphological quenching" process 
is truly important in maintaining the low observed rates of 
star formation in red sequence galaxies, we have performed 
the following experiment. We have binned galaxies with C < 
2.6 and C > 2.6 in the two-dimensional plane of NUV-r 
colour versus stellar mass (note that we use the same bin 
boundaries for both samples). We stack the Hi spectra of 
the galaxies in each bin and calculate the average Hi gas 
fraction, as explained above. In Figure [121 we report for each 
bin the ratio between the gas fraction of the disk-dominated 
(DD) objects and the bulge-dominated (BD) ones, i.e.: 




If the morphological quenching scenario is correct, then at 
fixed stellar mass, we would expect to find higher average 
Hi gas fractions for bulge-dominated galaxies on the red 
sequence than for disk-dominated galaxies on the red se- 
quence. 

The results in Figure [12] show that in general the op- 
posite is true. Gas fractions are always slightly higher for 
the disk-dominated galaxies than for bulge-dominated ones. 
The gas fraction differences do appear to be largest for red 
sequence galaxies with NUV-r > 5, but the sign o f the d if- 
ference contradicts the predictions of IMartig et a l. (2009). 

It is important to check that our result is not simply 
due to extinction effects. Some of the reddest, gas-rich disk- 
dominated objects may actually be heavily obscured systems 
that will move blueward when dust corrections are applied. 
Following Schi minovich et ahl (|2010T ) we have applied dust 
corrections to the NUV-r colours of the galaxies in our sam- 
ple with D n (4000)<l.fl In order to make up for the loss 
of red late- type galaxies, we had to apply stronger cut in 
concentration index (C>3) to define the early-type sample. 

5 D n (4000) is defined as the ratio of the average flux density in 
the continuum bands 3850 -3950 and 4000-4100 A, and traces the 
age of stellar populations jKauffmann et al.|[2003ah . 



The ratio of gas fraction between disk-dominated and bulge- 
dominated objects decreases, but we still find that the Hi 
gas fraction of the disk-dominated galaxies never falls be- 
low that of the bulge-dominated ones. In particular, if we 
divide the red sequence objects with 4.5 < NUV— r < 7 in 
two bins of stellar mass, we find a value of r=1.03 for the 
less massive galaxies, and r=1.71 for the more massive ones. 

We note that these results are in ag reement with 
those presented in ISchiminovich et al.l l|2010l ). In their pa- 
per, Schiminovich et al. use a volume-limited sample of 200 
galaxies from the GASS survey to explore the global scal- 
ing relations associated with the ratio SFR/Mhi, which they 
call the Hl-based star formation efficiency. They found that 
the average value of this star formation efficiency has little 
variation with any galaxy parameter, including the concen- 
tration index. 



6 SUMMARY 

We have carried out a stacking analysis using ALFALFA 
scans of a volume-limited sample of ~ 5000 galaxies with 
imaging and spectroscopic data from GALEX and the Sloan 
Digital Sky Survey. The galaxies have stellar masses greater 
than 10 10 M Q and redshifts in the range 0.025 < z < 0.05. 
We extract a sub-sample of 1833 "early-type" galaxies with 
inclinations less than 70°, with concentration indices C > 
2.6, and with light profiles that are well fit by a De Vau- 
couleurs model. We then stack the ALFALFA spectra of 
the galaxies from these two samples in bins of stellar mass, 
stellar surface mass density, central velocity dispersion, and 
NUV-r colour, and we use the stacked spectra to estimate 
the average Hi gas fractions Mhi/A/* of the galaxies in each 
bin. 

Our main result is that the Hi gas fractions of both 
early-type and late-type galaxies correlate primarily with 
NUV-r colour and stellar mass surface density. The relation 
between average Hi gas fraction and these two parameters 
is independent of C, and hence of the bulge-to-disk ratio of 
the galaxy. We note that at fixed stellar mass, early-type 
galaxies do have lower average Hi fractions than late-type 
galaxies, but this effect does not arise as a direct conse- 
quence of the presence of the bulge and we discuss possible 
implications below. 

We have also tested whether the average Hi gas con- 
tent of bulge-dominated galaxies differs from that of late- 
type galaxies at fixed values of NUV— r and /i*. We find 
no evidence that red-sequence galaxies with a significant 
bulge component are less efficient at turning their avail- 
able gas reservoirs into stars. This result is in contradiction 
with the "m o rphol ogical quenching" scenario proposed by 
IMartig et afl (|2009h . 



7 DISCUSSION 

We now consider possible implications of this work. 

1. The Hi content of a galaxy is independent 
of its bulge-to-disk ratio. This can be understood if the 
following two conditions are satisfied: 1) The Hi gas in early- 
type galaxies is always associated with disks, or with mate- 
rial that is in the process of settling into disks. 2) The for- 
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Figure 12. Morphological quenching test. For each bin of NUV— r and M* the ratio between the gas fractions of the disk-dominated 
galaxies with C < 2.6 f(HI)]j]j and the bulge-dominated galaxies with C > 2.6 f(Hl)BD is reported. The panel shows the adopted binning 
with the values measured for each bin. The arrows indicate upper/lower limit if one of the stacked spectra yielded a non-detection. Gray 
dots show how the targets (LTGs+ETGs) distribute along the plane. 



mation of the galactic disk is decoupled from the formation 
of the bulge. One way that condition (2) could be satisfied 
is if the disk is formed from gas that accretes well after the 
bulge formation event, so that the amount and the config- 
uration of the accreted material is not in any way related 
to properties of the bulge, such as its mass or its velocity 
dispersion. 

We now attempt to establish whether this hypothesis is 
correct. At fixed total stellar mass, we found that galax- 
ies with larger C have lower gas fractions than galaxies 
with smaller C (Figure [9] upper panel). If the gas is as- 
sociated mainly with the disk, then the ratio Mm /M*-,disk 
should not depend on C. We have fit the IWeinmann et al.l 
(2009) relation between B/T and C. Our bi-linear fit result 
is C = 1.19 + 2.28 ■ B/T. For each galaxy in our sample, 
we compute M*^. In Figure 13, we plot Mhi/M,;^ as a 
function of logM* for ETG and for sample A galaxies. As 
can be seen, the difference in gas fraction between the two 
samples is greatly reduced when the HI mass is divided by 
the disk stellar mass. For reference, the red and blue lines 
on the plot show Mui/M^m, as in Figures [7] and [9] (upper 
panels) . 

We note that our result that the Hi fraction of early- 
type galaxies depends on the size of the galaxy, but not on 
its bulge properties, is opposite to w hat is found f o r bulge 
stellar populations. In a recent paper, iGraves et ail ((2009b) 
carried out an analysis of 16,000 nearby quiescent galax- 
ies with 3-arcsecond aperture spectra from the Sloan Digi- 
tal Sky Survey. Their paper demonstrates convincingly that 
mean stellar age, [Fe/H], [Mg/H], and [Mg/Fe] scale strongly 



with the velocity dispersion of the bulge, but there is no de- 
pendence on R e at a fixed value of a. We thus infer that the 
star formation history of the central bulge does not depend 
on the size of the galaxy. Our own results show that the 
present-day gas content of early-types does depend on size. 
This again argues for bulge and disk formation processes 
that are decoupled. 

2. Galaxies with significant bulge component are 
not less efficient in turning their available gas into 
stars. One possible explanation of this result is that the rate 
of gas consumption in galactic disks is primarily regulated 
by externally-driven rat her than internally-g e nerat ed insta- 
bilities. In recent work, IChakrabarti fc Blitz! (|2009r i present 
an analysis of the observed perturbations of the Hi disk of 
the Milky Way and infer the existence of a dark sub-halo 
that tidally interacted with the Milky Way disk. In addition 
to dark sub-halos, luminous satellite galaxies are observed 
to interact with galactic disks. Finally, the dark matter envi- 
ronment at the cent ers of present day ha los is neither static 
nor in equilibrium. iGao fc White! (|2006T ) used the Millen- 
nium Simulation to study asymmetries in dark matter halo 
cores. They found that 20 per cent of cluster haloes have 
density center separated from barycenter by more than 20 
percent of the virial radius, while only 7 percent of Milky 
Way haloes have such large asymmetries. Because early-type 
galaxies reside in more massive dark ma tter halos than late- 
type galaxies (|Mandelbaum et alj|20 06). their disks may be 
subject to considerably larger externally-driven perturba- 
tions. All these effects may counteract the stabilizing effect 
of the bulge. 
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Figure 13. Average disk gas fractions of sample A (red circles) 
and ETG sample galaxies (blue circles) as a function of total 
stellar mass Af*. The lines show the results obtained measuring 
the total gas fraction, as in Figures [7] and [9] (upper panels). Gray 
dots show galaxies with ALFALFA detections from sample A. 



Further tests of the proposed external origin of the gas 
in early-type galaxies will come from more detailed analysis 
of its spatial distribution and kinematics. By studying the 
relation between gas and stellar angular momentum in early- 
type galaxies, one can hope to gain further understanding of 
how the gas was accreted. Studies of this nature are planned 
as part of the next generation of integral-field spectroscopy 
studies of nearby early-type galaxies, e.g. the ATLAS3D sur- 
vey (http://www-astro.physics.ox.ac.uk/atlas3d/ ). 



APPENDIX A: COMPARISON OF TWO 
DIFFERENT STACKING PROCEDURES 

In this section we compare two different methods for 
deriving mean Hi gas fractions from stacked spectra. 

1) Stacking Hi fluxes. In the first method, we derive 
the average Hi gas fraction using the equation 



,Mm_ _ 2.356 x 10 5 Dl((z)) E^L & ■ Wi 
y M» — 



1 + (z) (M») 



(Al) 



where (z) is the mean redshift of the stacked galaxies and 
(M t ) is their mean stellar mass. The main limitation is that 
this method may produce biased results if galaxies in the 
bin span a significant range in redshift or stellar mass. By 
stacking mainly non-detected spectra, we do not know how 
much each galaxy contributes to the total signal. This ap- 
proach will work best if the bins are small. We note that our 
redshift range is compact enough so that: 



Dl((z}) ~ (D L {z)Y 



(Dl(z).} 



(A2) 



We choose to split each stellar mass bin into three sub-bins 
in redshift: [0.025;0.033],[0.033;0.042],[0.417;0.050]. For each 
of these sub-bins, we evaluate the mean stellar mass (M*), 
and split the sample into two further bins according to: 
M*sS {Mt). We stack the spectra in each of the 6 sub-bins, 
measure a flux and evaluate a gas fraction using both the 



average and median values of z and M*. Then we evaluate 
the final gas fraction by averaging the six values obtained 
(weighted by the number of objects co-added). 

2) Stacking "gas fractions". This approach is 
described in section 13.21 and is the one that is adopted 
throughout the main body of this paper. 

First, we checked that the rms of the co-added spec- 
tra decreases as 1/vJV; when method (2) is applied, i.e. by 
multiplying the flux by distance and stellar mass, we are 
in fact rescaling the noise in the spectra in a sensible way. 
Results are shown in Figure [14] (panel a) : red circles are ob- 
tained by stacking "gas fractions" , black circles are obtained 
by stacking fluxes. The black points have been multiplied 
by the average squared distance and divided by the aver- 
age M* in order to have the same units as the red points 
[mjy Mpc 2 M" 1 ]. The dashed line is the expected 1/-/N de- 
pendence, which is still recovered using method (2). 

We then compared the results of stacking methods (1) 
and (2). We divided all the galaxies in sample A that were 
detected by ALFALFA (detection codes 1+2) into five stellar 
mass bins. We compare results from the different stacking 
methods with the results obtained by averaging together the 
individual catalogued measurements. The results are shown 
in Figure IT4l (panel b), where red circles show results ob- 
tained using method (2), and green triangles/blue squares 
are from method (1), when mean or median values of z and 
Mt. are used in equation (Al), respectively. The stars show 
the results from averaging together the individual detec- 
tions. Gray dots show the ALFALFA detections. The sig- 
nal recovered from the stacking is consistent with the mean 
value of the individual detections for each bin. The two dif- 
ferent stacking methods yield results that are also consis- 
tent with each other. Small differences of around ~ 10 % in 
Mhi/M* do occur in the two largest mass bins that contain 
the fewest objects. 

We also compared results for the different methods includ- 
ing the non-detections. In Figure [15] we show the same cor- 
relations studied in section U for sample A, computed using 
different methods. Symbols are the same as described above. 
Once again, we obtain good consistency. 
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Table Al. Average gas fractions for the samples shown in Figures 171 101 sample A (Fig. [7] and [8] red circles), the main ETG sample 
(Fig. [9] and [TO] blue circles), and the two ETG sub-samples defined with more stringent cuts on C and inclination (Fig. I9land ll0l black 
and cyan circles). For sample A we also report the slopes a (third column) of the relations log ( Mhi/M*} = a{ x) + b, where x is the 
quantity listed in the 1st column. 
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